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Abstract 

We derive mass corrections for semi-inclusive deep inelastic scattering of leptons from nucleons 
using a collinear factorization framework which incorporates the initial state mass of the target 
nucleon and the final state mass of the produced hadron. The formalism is constructed specifically 
to ensure that physical kinematic thresholds for the semi- inclusive process are explicitly respected. 
A systematic study of the kinematic dependencies of the mass corrections to semi-inclusive cross 
sections reveals that these are even larger than for inclusive structure functions, especially at very 
small and very large hadron momentum fractions. The hadron mass corrections compete with 
the experimental uncertainties at kinematics typical of current facilities, and will be important to 
efforts at extracting parton distributions or fragmentation functions from semi-inclusive processes 
at intermediate energies. 
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I. INTRODUCTION 



In recent years semi-inclusive deep inelastic scattering (SIDIS) has received much atten- 
tion as a tool to investigate various aspects of hadron structure, such as the flavor dependence 
of the nucleon's parton distribution functions, both unpolarized and polarized, through fla- 
vor tagging of hadrons in the final state. Observation of the momentum distribution of 
produced hadrons also allows access to the largely unexplored transverse momentum depen- 
dent parton distributions, which reveal a much richer landscape of the spin and momentum 
distribution of quarks in the nucleon, and which are the subject of increasingly greater focus 
at modern facilities such as Jefferson Lab. 

At high energies the scattering and hadronization components of the SIDIS process fac- 
torize and the cross section can be represented as a product of parton distribution and 
fragmentation functions. In practice, however, experiments are often carried out at few- 
GeV energies with Q 2 as low as 1 GeV 2 , suggesting that 1/Q 2 power corrections must be 
controlled in order to determine the applicability of partonic analyses of the data. 

One of the standard finite-Q 2 corrections that must be applied in analyses of inclusive 
deep inelastic scattering (DIS) data is target mass corrections (TMCs) Kinematical in 
origin, TMCs arise from leading twist operators in QCD, but enter as 1/Q 2 corrections to 
structure functions [2j. They are especially egregious at high values of the Bjorken scaling 
variable xb, even at relatively large Q 2 , and are crucial for reliable extractions of parton 
distributions in this region. To date, however, the phenomenology of TMCs has not been 
systematically considered in SIDIS, and we do so in this paper. 

Target mass corrections in inclusive DIS have usually been formulated within the operator 
product expansion, in which the subleading 1/Q 2 corrections arise from twist-two operators 
involving derivative insertions into quark bilinears ^, 0, 0, [f], Q, [s, 0]. Unfortunately, this 
method cannot be rigorously extended to the production of hadrons in the final state. An 
alternative approach to computing TMCs makes use of the collinear factorization (CF) 
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framework 



0, 
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12| , which has recently been used in both unpolarized |13|, [14|, ll5|] and 



polarized 16| inclusive DIS. Because here one works directly in momentum space, the method 
can be readily extended to SIDIS. In contrast to inclusive DIS, where the only mass scale 
entering the problem is that of the target hadron, in SIDIS finite-Q 2 corrections arise from 



both the target mass and the mass of the produced hadron. For generality we shall refer to 
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their combined effects as "hadron mass corrections" (HMCs). 

Hadron mass corrections in SIDIS at finite-Q 2 kinematics in CF were considered previ- 
ously in Refs. 



17 



in different collinear frames. Albino et al. [17j studied the effects of 



;he final state hadron mass, but did not consider the effects of the target mass. Mulders 
JjI derived corrections due to both target and produced hadron mass, but did not discuss 
the phenomenological consequences. Neither of these, however, addressed problems related 
to kinematic thresholds. 

In this work we use the CF framework to derive the mass corrections to the SIDIS 
cross section at finite Q 2 , and systematically investigate their implications at kinematics 
relevant to current experiments. The formalism is constructed specifically to ensure that 
physical kinematic thresholds for the semi-inclusive process are explicitly respected. In 
Sec. [Til we review the collinear formalism and discuss its application to semi- inclusive hadron 
production. To expose the origin of the corrections we work at leading order in a s ; next- 
to-leading order effects can be included in subsequent analyses. In Sec. [TTTJ we explore the 
relative importance of the HMCs numerically, and evaluate the size of the corrections in the 
cross sections and fragmentation functions at various kinematics. To assess their possible 
impact on data analyses, we also compare the magnitude of the HMCs at kinematics typical 
of modern facilities, such as Jefferson Lab and HERMES, with experimental errors from 
recent experiments. Finally, in Sec. IIVI we summarize our results and outline avenues for 
future developments of this work. A discussion of the formulation of HMCs in different 
collinear frames is presented in Appendix [A] 



II. SEMI-INCLUSIVE SCATTERING AT FINITE Q 2 

We begin the discussion of SIDIS at finite values of the photon virtuality Q 2 by defin- 
ing the relevant kinematics and momentum variables in a collinear frame, and introduce 
the hadronic tensor computed in a covariant parton model. Collinear factorization is then 
performed in the leading order approximation in which the produced hadron is effectively 
collinear with the scattered parton, which more directly reveals the effects of hadron masses 
on the cross sections and fragmentation functions. 
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A. External kinematics 



The four-momenta of the target nucleoli (p), virtual photon (q) and produced hadron h 
(ph) can be decomposed in terms of light-cone unit vectors n and n as 10] 

M 2 

p v = p + n v _j n v n & \ 

2p + 

(f = -iP + n^ + ^-n^ (lb) 

u £, m h± H u , ChQ u a , n. 

where M is the target nucleon mass, Q 2 = —q 2 , and the light-cone vectors satisfy n 2 = n 2 = 
and n -n — 1. Here we define light-cone components of any four- vector v by v + = v ■ n = 
(v o + f z) / v2 and v~ = v ■ n = (vq — v z ) / y/2. The momenta p and q are chosen to lie in the 
same plane as n and n, as for inclusive DIS. We call this the (p, q) collinear frame; other 
possible choices are discussed and compared in Appendix A. The nucleon plus-momentum 
p + can be interpreted as a parameter for boosts along the z-axis, connecting the target 
rest frame to the infinite-momentum frame; the target rest frame (p + = MJ y/2) and the 
Breit frame (p + = Q/{y2£)) are part of this family of frames. The transverse momentum 
four- vector of the produced hadron p% ± satisfies ph± ■ n = ph± ■ n = 0, and we define the 
transverse mass squared as = — p1 ± , where rrih is the produced hadron mass, and 
the transverse four- vector squared is p 2 h± = —pl±- 

In the chosen collinear frame the variable £ — —q + /p + defined in Eq. fllbl) coincides with 
the finite-Q 2 Nachtmann scaling variable 0, [3] , 

£ = 2xb ( 2 ) 

l + y/l + 4x 2 B M 2 /Q 2 ' 

which in the Bjorken limit (Q 2 — > oo at fixed xb) reduces to the Bjorken scaling variable 
%b = Q 2 /2p-q. The scaling fragmentation variable (h = Vhl^T defined in Eq. (fTc|) is related 
to the fragmentation invariant Zh = Ph • p/q • P by 




Ax 2 B M 2 m 2 h± 



(3a) 



A Q 4 

and the positivity of the argument in the radical in Eq. ( |3al) is ensured by the condition 
Eh > rrihx, which imposes 

z h > zr = 2x B ^ . (3b) 
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One can also define (h in terms of the invariant T]h = 2ph ■ q/q 2 by 

Vh 

2 V" V" " 



a=^(l + A/H-3^) , (3c) 



which is convenient for discriminating between the target and current fragmentation hemi- 
spheres in hadron production. Note that in the target rest frame Zh — E^/v is the usual 
ratio of the produced hadron to virtual photon energies. In the Breit frame 77^ = Ph z /Qz is 
the ratio of the longitudinal components of the hadron and photon energies, which can be 
used to define the current (r] h > 0) and target (i] h < 0) hemispheres for hadron production. 
In the Bjorken limit one has (h Zh — > f]h- 

Conservation of four- momentum and baryon number impose an upper limit on the Xb 
variable, 

*B<(l + m;+ Jf** )'W, (4) 

which corresponds to the exclusive production of a nucleon and a hadron h in the final state. 
Similarly the limits on the fragmentation variable (h are given by 

<Ch<l + £7yr, (5) 



1 _ £ Q2 - ™ - - Q2 

where the lower limit corresponds to diffractive production of the hadron h, and the upper 
limit reflects the fragmentation threshold, which approaches unity in the Bjorken limit. 

B. Parton kinematics in collinear factorization 

At the partonic level the SIDIS process at leading order in the strong coupling constant 
a s is illustrated in Fig. [H It proceeds through the scattering from a quark carrying a 
light-cone momentum fraction x = k + /p + , which then fragments to a hadron h carrying a 
light-cone momentum fraction z = p^/l" , where k and I are the four- momenta of the initial 
and scattered quarks. At higher orders the hard scattering can also take place from a gluon, 
and additional partons can be created in the collision. 

The parton momenta k and I can be parametrized in terms of the light-cone vectors n 
and n as 

k 2 4- k 2 

k" = xp + + x + h!\ , (6a) 

2xp + 

72 4. 72 „- 

l» = -—^ n» + ^n» + l» (6b) 

2p h jz z 




FIG. 1: Kinematics of semi- inclusive deep inelastic lepton-nucleon scattering at leading order, 
producing a final state hadron h. The momenta of the target nucleon (p), virtual photon (q), 
incident (k) and scattered quarks (/), and the produced hadron (ph) are labeled explicitly, with 
and A denoting the correlators relevant to the quark distribution and fragmentation functions. 
The vertical dashed line represents the cut of the forward amplitude. 

with the parton transverse momentum four-vectors k± and l± orthogonal to n and n. 
In collinear factorization the hard scattering amplitude is expanded around on-shell and 
collinear momenta k and I, 

k 2 

P = X p + n^-\ -n M (7a) 

2xp + 

p = P+PhJz 2 ^ + P±^ + PL ( 7b ) 
2 Ph /z z z 

where the initial and final collinear parton "masses" k 2 and I 2 are kept for generality. 

Denning the invariant x = —q 2 /2k ■ q as the partonic analog of the Bjorken variable xb, 
at finite Q 2 one has 



Using the methods described in Ref. 15|] one can show that for SIDIS cross sections inte- 
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grated over ph±, x is constrained to be in the range 

m\ k 2 f H ^m 2 h \ ^ 1 1 ( _ 2Mm h + k 2 



where the lower limit arises from the minimum of the current jet mass, and the upper limit 
corresponds to collinear spectators with minimal mass. These limits agree with the limit on 
xb in Eq. (j4]) for any k > x((h — 1)Q 2 /£• For the fragmentation process one finds analogous 
limits on (h, 



a <k* <i + !^> (io) 



which agrees with the limit in Eq. (jSJ), provided that k 2 < xM 2 . The requirement that 
the collinear parton masses be independent of the parton momentum (viz., independent of 
x) implies k 2 < 0. Combined with the above lower limit on k 2 , this naturally leads to a 
collinear expansion around a massless initial state parton, k 2 = 0. 

The choice of I 2 is made by considering the cross section at leading order in a s . Four- 
momentum conservation for the hard scattering, together with the choice k 2 = 0, leads to 
the relations x — £(1 + l 2 /Q 2 ) = £h an d z = Ch- Clearly z falls within the kinematic limits 
ffTUj) . However, in order for x to respect the limits ([9]) we choose I 2 = m\j^ in which case 

While larger values of I 2 would also allow x to fall within the limits ([9]), this choice is the 
closest to the physical quark mass. 

We stress that our prescription for the collinear parton masses k 2 and I 2 is dictated by the 



external kinematic limits in Eqs. (j3J) and (jSJ), which are independent o 
and collinear factorization approximations. As discussed in Refs. 15 



the parton model 



21], this is crucial 



when considering cross sections close to the kinematic limits, such as at large xb or large 
Zh- However, as we shall see in the next section, the SIDIS cross section can also receive 
non- negligible corrections at small xb since ^ > £ ~ xb- This is qualitatively different from 
the behavior of the target mass corrections in inclusive DIS, which are always suppressed at 
small xb 
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C. Hadron tensor and cross section at leading order 



In collinear factorization the hadron tensor at leading order, to which we restrict the rest 
of our analysis, can be written as 

2MW^(p,q,p h ) = [ d A k d A l 5 {A \k + q-l) Tr[%(p, k) ^ A h q (l,p h ) 7 "] , (12) 

where the sum is taken over quark flavors q, and the correlators and A^ e ncode the rele- 
vant quark distribution and fragmentation functions, respectively [ll|, H, 18]. According to 
our prescription for the collinear momenta, the 5-function depends on the collinear momenta 
k and I, so that integrations over dk~ d 2 k± and dl + d 2 l± act directly on the correlators $ 
and A. The leading twist part of the cross section can then be extracted by retaining the f£ 
and ifi components in the Dirac structure expansion of the integrated correlators, 

J dk~d 2 k ± $,(p, k) = ~f q {x)$ + ... , (13a) 
J dl + d 2 l ± A h q (l,p h ) = \D\{z)i + . . . , (13b) 



where the dots indicate contributions of higher twist [20j . The nonperturbative quark dis- 
tribution function f q (x) and quark-to-hadron fragmentation function D q (z) are explicitly 
defined as 



/«,(*) = \j dk-d 2 k L Tr [l + %(p,k Jlk+=xp 
lc 1 f dw 



e^"(iV|^(0)7 + ^(^-n)|iV) , (14a) 



2 ./ 2tt 

=P h I* 



D h q {z) = Z - J dl+d'hTrfr-A^p^ 

= f E/ ^e^/^ + (0|^( U ; + n)|/ i ,X)(/ i ,X|^(0)7"|0) , (14b) 



where "LC" denotes use of the light-cone gauge, and the fragmentation function is normal- 
ized such that J2 h Jg 1 dz z DHz) = 1 fis| . 

From Eq. (1121) the energy-momentum conserving 5-function can be decomposed along 
the plus, minus, and transverse components of the light-cone momentum. The plus and 
minus components yield a product of ^-functions that fix x = £/, and z = ( h , while the 
transverse component constrains the transverse momentum of the scattered quark to vanish, 
which restricts the produced hadrons to be purely longitudinal, ph± = z l± = 0. Hadrons 
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with nonzero transverse momentum can be generated from higher order perturbative QCD 
processes, or from intrinsic transverse momentum in the parton distribution functions, as 
in the case of transverse momentum dependent distributions [2Cj , but are not considered in 
this work. The resulting hadron tensor in the presence of hadron mass effects, 



2MW^(p,q,p h ) = j^ej S^(p ± )Tr[$^iY}f q ^ h )D h q (C h ), (15) 

<? 

is then factorized into a product of parton distribution and fragmentation functions evaluated 
at the fmite-Q 2 scaling variables ^ and (h, instead of xb and as would be obtained in the 
massless case, and recovered from Eq. (fTBT) in the Bjorken limit. Note that this prescription 
is the same as that used in Ref. 15] when discussing inclusive DIS in the presence of jet 
mass corrections, and is close in spirit to that advocated in Ref. 21[, where the trace is 
calculated as in the massless case, but overall parton momentum conservation respects the 
external kinematics. 

Finally, the SIDIS cross section is computed by contracting the hadron tensor with an 
analogous lepton tensor [2o|, leading to 

_ da 2na 2 y 2 d( h ^ 2 2 h 2 

° = dx B dQ 2 dz h = -QTT^-edz;^ f ^ Q)D ^ h ' Qh (16) 



where the dependence of the functions on the scale Q 2 is made explicit, and the Jacobian 
dC h /dz h = (1 - M 2 £ 2 /Q 2 )/(l - £, 2 M 2 m 2 /( 2 Q 4 ). In Eq. <^> the variable y defined as 
y = p ■ q/p ■ Pi, where pe is the lepton momentum, represents the fractional energy transfer 
from the lepton to the hadron in the target rest frame (y = v/E, with E the lepton energy), 
and e — (1 — y — y 2/ ~y 2 / 4) / (I — y + ?/ 2 [l/2 + 7 2 /4]) is the ratio of longitudinal to transverse 
photon flux, with 7 2 = Ax 2 B M 2 /Q 2 . The cross section differential in r/^ can be obtained 
using d(h/di]h = 1/(1 + m^/C^Q 2 ) instead of dC^h/dzh- It is interesting to observe that since 
£/j depends explicitly on nih and (h depends on Zh and xb, at finite Q 2 the scattering and 
fragmentation parts of the cross section ( TL6{) are not independent. 

As a final remark we note that at the maximum allowed xb for SIDIS, Eq. (j4j), the value 
of £/i is smaller than ^(xb = x^ ax ) < 1. As in the case of inclusive DIS [la ], the SIDIS cross 
section therefore does not vanish as xb —> x^ ax , which is a manifestation of the well-known 
threshold problem 1|. On the other hand, from Eq. (TTUT) the fragmentation variable (f^ < 1, 
and no threshold problem appears in the fragmentation function since D((h) — > as ^ — > 1. 
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FIG. 2: Ratio of cross sections a/a^ for semi-inclusive charged-pion production ((vr + + 7r )/2) as 
a function of at several Q 2 values for (a) xb = 0.3 and (b) xb = 0.8. 

In the next section we shall examine the phenomenological consequences of the fmite-Q 2 
rescaling of the SIDIS cross section numerically. 



III. HADRON MASS CORRECTIONS 



Using the hadron mass corrected expressions for the SIDIS cross section derived above, 
we next explore the dependence of the cross sections and fragmentation functions on the 
fragmentation variable Zh, for various xb and Q 2 values and for different final state hadron 
masses. We then compare the relative size of the HMCs with the experimental uncertainties 
from recent SIDIS experiments at Jefferson Lab and the HERMES Collaboration, as well as 
with higher energy data from the European Muon Collaboration (EMC) and HERA. 



A. HMC phenomenology 

To illustrate most directly the effects of the HMCs, in Fig. [2] we consider charged pion 
production (average of 7r + and n~) and plot as a function of z h , for different xb and Q 2 , the 
ratio of the full cross section a in Eq. (jT6l) to the cross section a^°\ defined by taking the 
massless limit for the scaling variables 

(£h —> x B , Ch — > z h ) and setting d( h /dz h = 1 



22 



For the numerical computations we use the leading order CTEQ6L parton distributions 
and the KKP leading order fragmentation functions [23[], unless otherwise specified. The 
ratio at xb = 0.3 in Fig. [2^ a) is enhanced by < 20% at Q 2 = 2 GeV 2 for z h < 0.7, but rises 
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FIG. 3: Ratio of the hadron mass corrected isoscalar parton distribution function f q (£h) f° r 
q = u + d to the massless limit distribution fq^ as a function of xb, for irih = m n and Ch = 0.2. 



dramatically as z^ — > 1. The effect is naturally smaller at higher Q 2 values, but the rise at 
high Zh is a common feature for all kinematics. The same ratios at xb = 0.8 in Fig. EJ^b) 
show approximately an order of magnitude larger overall effect (note the logarithmic scale!). 

The small upturn in the ratios at low Zh for the lowest Q 2 in Fig. [2] can be understood from 
the interplay between the finite-Q 2 kinematics and the shape of the fragmentation function. 
Assuming the fragmentation function is smooth, one can expand the ratio of corrected to 
uncorrected functions in a Taylor series as 

D(z^~ + D(z^ {Ch ~ hh (17) 
The Zh dependence of the HMCs arising in the fragmentation function is mostly determined 
by the negative shift in the fragmentation variable ((h — Zh) an d by the local rate of change 
over Zh of the fragmentation function. The pion fragmentation function generally behaves 
as a negative power of Zh at small Zh, and the negative slope drives the ratio of corrected to 
uncorrected fragmentation functions upward as Zh — » ^™ m , where \(h — Zh\ is maximum. For 
kaons and protons the slope of the form factor can be positive, which would suppress the 
mass corrected cross section in the vicinity of £™ m . In the limit — *■ 1, on the other hand, 
the ratio a/a^ becomes divergent for any kinematics and any hadron species because the 
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FIG. 4: (a) Dependence of the ratio of SIDIS cross sections a/a^ with and without HMCs for 
different produced hadrons, h = ir + + tt~ , K + + K~ or p + p. (b) Ratio of cross sections for 
h = tt + + 7r~ for different values of the pion mass, relative to the massless cross section. In both 
cases the kinematics chosen are xb = 0.3 and Q 2 = 5 GeV 2 . 

cross section oc D(zh) vanishes, while the rescaled cross section remains finite. 

At very small values of Zh the factor (1 + Tn^/QhQ 2 ) in the definition of ^ in Eq. ffTTl) can 
render ^ h larger than xb, suppressing the ^-rescaled parton distributions relative to their 
asymptotic limit and driving a/a^ slightly below unity. As discussed below, for heavier 
hadrons this effect will be more pronounced. The effect of the £/, rescaling on the SIDIS 
cross section is illustrated explicitly in Fig. [3j where we show the ratio of the isoscalar parton 
distribution functions f q ,q = u + d, with [f q = f q (£h)] an d without [fq ^ = / 5 (xb)] hadron 
mass corrections, as a function of xb for Ch = 0.2 and rrih = fn n . At Q 2 = 2 GeV 2 the mass 
corrected parton distribution is several times larger than the uncorrected one at xb = 0.8, 
and even at Q 2 = 20 GeV 2 the HMC is some 50%, with the effect increasing dramatically as 
Xb — * I- This sharp rise is analogous to that in inclusive DIS 151 ] . and arises from ^ being 
smaller than xb when the latter is large. This is responsible for the large overall magnitude 
of the corrections in Fig. EJ^b) compared with those at xb = 0.3. In contrast, the ^ rescaling 
effect becomes quite small at xb J$ 0.3 for all the Q 2 considered, and in fact drives the ratio 
below unity, as discussed above. 

The relative importance of HMCs for different produced hadron species is illustrated in 
Fig. H](a), where the ratio a/a^ is shown as a function of Zh for xb = 0.3 and Q 2 = 5 GeV 2 . 
Over the range 0.3 < z h < 0.8 the HMCs yield an upward correction of < 10% for the pions, 
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FIG. 5: (a) Ratio of hadron mass corrected to uncorrected fragmentation functions D/D^ for 
favored (solid) and unfavored (dotted) production of 7r + , for = 0.8 and Q 2 = 2 GeV 2 . (b) 
Comparison of the hadron mass correction 5 HMC i? to the ratio of unfavored to favored fragmen- 
tation functions R = D~ / D + with experimental errors 5 cxp R on R from the recent Jefferson Lab 
experiment E00-108 [25], normalized to the central values of the data points. 



but a downward correction of < 20% and < 40% for kaons and protons/antiprotons, respec- 
tively. At lower Zh the cross section ratio for the heavier hadrons decreases dramatically 
because of the large suppression of the parton distribution from the (1 + m 2 h /(,hQ 2 ) factor 
in which overwhelms any other small- ^ effect. 

Note that in Fig. H^a) the appropriate fragmentation function for each produced hadron 
species has been used, which introduces a flavor dependence in the HMC because of the 
different fragmentation function shapes for each hadron. To isolate the effects of the hadron 
mass alone, in Fig. H](b) the ratios of cross sections computed with charged pion fragmenta- 
tion functions and masses = m n (= 0.139 GeV), 0.5 GeV and 1 GeV are shown relative 
to the cross section with = 0, for which (h = Zh^/xs- One can see that in general 
increasing the hadron mass suppresses the cross section because of the £/, scaling, and the 
inversion of the HMC hierarchy in Fig. Hl(a) going from low to high Zh is due to the increas- 
ingly negative large- z slope of the fragmentation functions for kaons and protons. While the 
differences at the physical pion mass are very small, for larger hadron masses ~ 1 GeV the 
effects can be quite significant at < 0.4 even for Q 2 values of several GeV 2 . 



13 



B. Implications for experiments 



One of the unique capabilities of SIDIS is the ability to tag individual quark flavors by 
selecting specific hadrons in the final state. For example, because of its valence quark con- 
tent, production of 7r + is mostly sensitive to the u quark, requiring only a single qq pair 
creation from the vacuum, while 7r~ reflects mostly the d quark content of the target nu- 
cleon. This simple picture of primary fragmentation provides a good approximation to the 
production mechanism only at large Zh, however, and at low Zh secondary fragmentation 
involving two or more qq pair productions dilutes the direct flavor tagging. The primary 
fragmentation process is parametrized by the "favored" fragmentation function D + , describ- 
ing u — ► 7r^ or d — > 7i~ hadronization, while the secondary fragmentation is parametrized by 
the "unfavored" fragmentation function D~ , describing u —>■ n~ or d —>■ n + hadronization. 

Because the D + and D~ functions have rather different z^ dependence, with unfavored 
fragmentation strongly suppressed at large Zh, they will be affected differently by the hadron 
mass corrections: one would expect larger HMCs for the unfavored process since the mag- 
nitude of the slope \dD(zh)/dZh\ in Eq. (117)) is larger for D~ than for D + . In Fig. 0a) one 
observes precisely this; here, we provide an upper limit (given the choice of xb = 0.8 and 
Q 2 = 2 GeV 2 ) to the relative size of the mass effect in D/D^°\ which is universally larger 
in the unfavored fragmentation function. In the numerical computations we have used the 
favored and unfavored fragmentation functions from Ref. 24j. At lower xb the correction 
will be smaller, although the qualitative features of the effect will remain. 

The relevance of the HMCs to experimental data on the ratio R = D~ / D + is expressed 
in Fig. 0(b), which directly compares the difference 5 nMC R = (D~ / D + ) — (D~ / D + )^ to 
the statistical uncertainty S exp R in the extracted values of R from the recent Jefferson Lab 
experiment E00-108 25| at xb = 0.32 and Q 2 ~ 2.5 GeV 2 , with both quantities normalized 
to the central values of the measured R ratio. While at small z^ the HMC is relatively 
small compared with the experimental errors, at large Zh (> 0.6) it begins to compete with 
the experimental uncertainty, suggesting that the hadron mass here poses a non-negligible 
effect. 

The importance of the hadron mass corrections for experimental cross sections is ex- 
amined in Fig. [61 where we compare the calculated difference 5 nMC a = a — with the 
experimental uncertainties 5 exp a, normalized to the central values of the cross section for 
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FIG. 6: Comparison of the hadron mass correction to the SIDIS cross section for charged hadron 
production, <5 HMC cr, relative to the experimental cross section, S cxp a, with the relative experimental 



26J] at Q 2 = 2.5 GeV 2 and x B = 0.082, 



uncertainty as a function of Zh for (a) HERMES experiment 
and (b) Jefferson Lab experiment E00-108 [25] at a similar Q 2 but at xb = 0.32 



charged hadron production from HERMES 



261 ] and Jefferson Lab 25j. Since both exper 



iments are dominated by the semi-inclusive production of pions, so that £h ~ £, HMCs 
generally produce upward shifts relative to data. The mass corrections at the HERMES 
kinematics in Fig. E(a), where Q 2 ~ 2.5 GeV 2 and xb = 0.082, are generally very small 
compared with the size of the experimental uncertainties. At higher energies, HMCs to 
fixed-angle measurements by the EMC 27| at large xb values are also found to be negligible 
due to suppression by Q 2 , which increases with xg. Were these experiments conducted at 
smaller angles, however, it is likely that HMCs would become important. 

On the other hand, for the Jefferson Lab experiment E00-108 [^J in Fig. EJ^b) at a 
similar Q 2 but larger xb = 0.32, the mass effects are approximately 2 times larger than 
the experimental statistical errors. This illustrates the potentially significant impact that 
HMCs can have on leading-twist analyses of SIDIS data at moderate and large xb and low 
Q 2 . To avoid these effects one would either need to go to smaller xb or larger Q 2 values, for 
example afforded by the 12 GeV energy upgrade at Jefferson Lab. Alternatively, since the 
HMCs are calculated and model independent, lower Q 2 and higher xb data will still yield 
useful leading twist information provided the mass corrections are accounted for. 

Measurements at small xb ~ 0.001 and Q 2 > 12 GeV 2 have been performed by the HI 
collaboration 31| at HERA, and the data presented in terms of the fragmentation invariant 
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r]h- The phenomenology of HMCs is markedly different in terms of r/h from that discussed 
thus far in terms of Zh because of the different functional forms for (h in Eqs. ([3]), which 
constrains (h > Vh, and because of the Jacobian d(h/dVh- In their analysis of the HI data, 



Albino et al. 



171 | included the (h rescaling of the fragmentation process, but neglected the 



effects of the target mass, which would be problematic for heavier hadrons such as kaons 
and protons. The HI Collaboration measured charged hadron multiplicities, dominated 
by pions (~60%), with smaller contributions from kaons (~30%) and protons (~10%). In 
the measured Q 2 range the m^/Q 2 term in £/, is therefore strongly suppressed and at the 
typically low x_r values one has £ ~ xb, so that overall we find the HMCs to be similar to 



those in Ref. [l7j. However, for identified kaons, and especially protons, the SIDIS cross 
section would be more strongly suppressed compared to the results of Ref. [17] because 
£h ~ xb{^ + m h/Q 2 ) is significantly larger than xb- This suppression may be non-negligible 
for the extraction of kaon and proton fragmentation functions from small-Xe data. 



IV. CONCLUSIONS 



In this paper we have derived hadron mass corrections to semi-inclusive deep inelastic 
cross sections at finite Q 2 and have performed a systematic exploration of their phenomeno- 
logical consequences. Within the collinear factorization framework the modifications to the 
leading order SIDIS cross sections from initial and final state masses arise from a rescaling 
of the quark distribution and fragmentation functions in terms of the modified Nachtmann 
scaling variable ^ and a finite-Q 2 fragmentation variable (h, respectively. The need for a 
modified Nachtmann variable is dictated by the requirement that the physical kinematic 
thresholds for the semi- inclusive process are explicitly respected. 

We have examined the effects of the hadron mass corrections numerically at kinematics 
relevant to recent experiments, finding sizable effects at both small and large values of Zh, 
as well as for increasing xb and m^, and low Q 2 . Our results emphasize the importance 
of controlling for such corrections in intermediate to high-x^ experiments executed at low 
Q 2 , of which measurements at Jefferson Lab are typical, although not exclusive. We find 
that the hadron mass corrections can in some cases compete with the quoted experimental 
errors (as in the measurement of the ratio of unfavored to favored fragmentation functions 
D~/D + ) or overwhelm them (as in the cross section measurements). Due to the presence 
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of the modified Nachtmann variable the HMCs may also need to be considered for small-x^ 
collider experiments for the production of heavier hadrons such as kaons and protons, and 
their effects require further study. 

The most direct use of the results presented here will be in leading twist analyses of 
SIDIS cross sections, where the HMCs must be included before extracting information on 
parton distribution and fragmentation functions, especially at large xb and z^. Application 
of this work can also be found in studies of semi-inclusive data in the nucleon resonance 
region, which has been the focus of attention recently in view of understanding the workings 



of quark-hadron duality [28|, |29|, |30|] . 



While the present analysis has been performed at leading order in a s , in the future we 
plan to extend the formalism to next-to-leading order, which will permit a more quantitative 
treatment of transverse mass dependence of the produced hadrons, ph± 7^ 0. It will also allow 
contact with transverse momentum dependent parton distributions, in which nonzero parton 
transverse momentum, k± 7^ 0, is an essential element. Finally, as in the inclusive DIS case, 
the SIDIS cross section corrected for hadron mass effects exhibits the threshold problem 
which renders it nonzero as -> x^ x . Solutions of this problem proposed in the literature 
for inclusive structure functions 4], 0, [ijj] will be extended to SIDIS in future work. 
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APPENDIX A: COLLINEAR FRAMES IN SEMI-INCLUSIVE DIS 

A collinear frame in Minkowski space is defined by any two four- vectors. The intersection 
of the plane where they lie with the light-cone defines the light-cone four- vectors rT M and n^, 
that satisfy n 2 = n 2 = and n ■ n — 1. In SIDIS the hadronic tensor depends on the three 
vectors p M , g M and p 1 ^, which define three possible collinear frames denoted (p, q), (ph, q) and 



17 



iPhiP)- The (p,q) frame is the only frame that can be defined in DIS and is the one used 
in this work; the (ph, q) frame is the only one that can be defined in semi- inclusive hadron 
production in e + e~ collisions; and finally the (ph, p) frame is typically preferred for analysis 
of transverse momentum dependent parton distributions in SIDIS. 

In terms of the vectors p, q and ph one can define two fragmentation invariants, 

Ph ■ P 2p h -q 

Zh = , Vh = 5— , (Al) 

q ■ p q z 

which together with xb, Q 2 , M 2 and ml form a complete set of scalar Lorentz invariants in 
SIDIS. Because the variable r] h is defined independently of the target momentum, the effects 
of the final state hadron mass will decouple from those of the target mass in all reference 
frames. In contrast, is defined in terms of both the target and produced hadron momenta, 
so that the target and hadron mass effects here will be entangled. 

The light-cone fractional momentum £ (Nachtmann variable) and the fragmentation vari- 
able (h are defined in terms of the plus and minus components of the momenta as in Eqs. (T2J) 
and f[3"a|) . 

£ = — , Cfc = - ■ (A2) 

p + q~ 

We use these definitions in all three frames; however, in each frame the light-cone vectors 
(and therefore the plus and minus components of the four-momenta) will be different. In 
the following we discuss each of the three collinear frames and the consequences within each 
frame of the choice of fragmentation invariant. 

(p, q) frame. In this frame the external vectors can be decomposed in terms of the 
light-cone vectors n and n as in Eqs. (00), 

M 2 

p» = p+n 11 + — -n M , (A3a) 
2p + 

(f = -£ p +n» + ^f-n» , (A3b) 
m 2 O 2 

A = i^re + * +ih h n * +A ' (A3c) 

where m\ ± = rn 2 h —p 2 h± = Tn 2 h +p 2 hl _ is the transverse mass of the produced hadron. Inverting 
the definition of xb, the Nachtmann scaling variable can be written as in Eq. (121) , 

; ^ (A4) 

l + y/T+ 4x 2 B M 2 /Q 2 



18 



Similarly, the hadron fractional momentum (h can be expressed in terms of either the frag- 
mentation invariant Zh, 



z h £ ( i / 1 A x 2 b M 2 m 2 

U = 

or in terms of rjh 



^y^ + V 1 - 4 ?^ 1 ' (A5) 



2 



* = T ( 1 + V 1 + 4 ^' ' (A6) 

One can show that for any finite (h the variable z^ — > f]h in the Bjorken limit. This is 
obviously true in any frame. 



(Phi q) frame. In this frame, used in Ref. [17j for example, the external SIDIS vectors 
are defined as 

M 2 

p k = p +n» + -±71^ + pi , (A7a) 
zp 

(f = -fr+n» + -^-n" , (A7b) 
2£p+ 

A = ^ p+ * +i ^ n '' (A7c) 

where M]_ = M 2 — p\ = M 2 + p\ is the transverse mass of the target nucleon. The 
Nachtmann variable in this case is given by 

= 2x B 

l + y/l + AxlM'i/Q 2 ' 

which, in contrast to its definition in the (p, q) frame, depends explicitly on the transverse 
mass of the target nucleon. Furthermore, in terms of the fragmentation invariant Zh, the 
finite-Q 2 fragmentation variable (h is given by 



^H i - 4 ¥P ] • (A9) 



or in terms of rjh by 



C^fl + Jl+^f ] • (A10) 
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(jPhiP) frame. The external vectors in this frame, used in Ref. 18( for example, are 
given by 

M 2 

p» = p+n^ + , (Alia) 

Zp 

<f = -fr+n* + + q£ , (Allb) 

A = A^ +Ch &'- (Auc) 

where Q\ = Q 2 — q]_ = Q 2 + Q± is the transverse mass of the virtual photon. The Nachtmann 
variable in this frame depends explicitly on Q 2 ±} 

^=91 2x * (A12) 

and the fmite-Q 2 fragmentation variable is given by 



z h £ Q 2 / / x 2 M 2 m 2 i 



or 



Vh Q 2 L, L imlQ 2 

^m 1 + f + nl-Q 



Relations between frames. In general the frames discussed here are distinct. However, 
to leading order in 1/Q 2 the vectors p, q and ph lie in the same plane and the three frames in 
fact coincide. Comparing the (p, q) and (ph, q) frames, for example, the differences between 
the transverse momenta and scaling variables can be expressed as 

Ph± =Pl + 0(p* 2 /Q 2 ) , (A15a) 
£ =C + 0(p* 2 /Q 2 ), (Al5b) 

a =c h +o( P f/Q 2 ), (Aisc) 

where the asterisks (*) label quantities in the (ph, q) frame. Similar relations are applicable 
also for the parton fractional momentum x and the hadron fractional momentum z. At 
leading order in collinear factorization one has Ph± = and the frames are manifestly 
equivalent. Moreover, since (p 2 ± ) <C Q 2 for p/^-integrated cross sections at next-to-leading 
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order the differences between the collinear frames should remain small. It will nevertheless 
be important to check whether, and in what kinematic range, this approximation is valid. 
On the other hand, for unintegrated cross sections the differences between frames become 
relevant and their effects must be quantified. 
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